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Abstract: The many-body ground state of a very general class of electron- 
phonon Hamiltonians is proven to contain a spin singlet (for an even number of 
electrons on a finite lattice). The phonons interact with the electronic system 

in two different ways — there is an interaction with the local electronic charge 

O 



ctf 



and there is a functional dependence of the electronic hopping Hamiltonian on 
the phonon coordinates. The phonon potential energy may include anharmonic 
terms, and the electron-phonon couplings and the hopping matrix elements may 
be nonlinear functions of the phonon coordinates. If the hopping Hamiltonian 
is assumed to have no phonon coordinate dependence, then the ground state is 
also shown to be unique, implying that there are no ground-state level crossings, 
and that the ground-state energy is an analytic function of the parameters in 
the Hamiltonian. In particular, in a finite system any self-trapping transition is a 
smooth crossover not accompanied by a nonanalytical change in the ground state. 
The spin-singlet theorem applies to the Su-Schrieffer-Heeger model and both the 
spin-singlet and uniqueness theorems apply to the Holstein and attractive Hub- 
bard models as special cases. These results hold in all dimensions — even on a 
general graph without periodic lattice structure. 
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I. INTRODUCTION 

Electrons have a tendency to pair when the effective electron-electron interaction has 
an attractive region; in particular this occurs when electrons interact by exchanging bosons. 
The resulting ground state then often sustains either superconducting or charge-density- 
wave order. The simplest interacting Hamiltonian of this type is one in which electrons 
interact indirectly with each other via phonons. Migdal 1 analyzed the electron-phonon in- 
teraction in the normal state and discovered that in the limit in which the phonon frequency 
O is much smaller than the Fermi energy Ef, the full many-body theory can be described 
by a first-order, self-consistent Hartree-Fock theory, and that the neglected higher-order di- 
agrams (vertex corrections) usually contribute to order Q/Ef. This result has been named 
Migdal's theorem and it classifies those nonadiabatic processes that are typically important 
for describing low-frequency electron-phonon interactions. Soon thereafter, Eliashberg 2 
generalized Migdal's result to the superconducting phase and discovered that a similar 
first-order self-consistent Hartree-Fock theory would describe superconductivity. Rowell 
and McMillan 3 subsequently demonstrated that one could directly measure the electron- 
phonon spectral density from tunneling experiments and then use the formalism of Migdal 
and Eliashberg to describe all of the remaining properties of the superconducting state. 
Migdal-Eliashberg theory has been successful in predicting transition temperatures (and 
other materials properties) of most low temperature superconductors 4 ' 5 . 

The electron-phonon Hamiltonian considered here is 

o x,yEA aieA j = l •> 

(1.1) 

We can also add to this an attractive Hubbard type interaction --as discussed later in 
Sect. III. Our notation is the following: The electrons occupy positions on a finite "lattice" 
or "graph" A, which is some collection of |A| sites; we emphasize that no specific periodicity 
or dimensionality is assumed. The operator c xa (c xa ) is a creation (annihilation) operator 
for an electron at lattice site x with ^-component of spin a =| or |. These operators satisfy 
the anticommutation relations c xa c ya + c ya c xa = 5 xy and c xa c ya + c ya c xa = for each 
a =| or [. It is customary to assume that the up-spin operators also anticommute with 
the down-spin operators, but it is more convenient for us to assume that they commute, 
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i.e., c x ^Cyi — Cyic x ^ = 0, etc.. This change is innocuous [as long as particle number is 
conserved -- which it is with the Hamiltonian (1.1)] and is effected by replacing c x -\ by 



1 r • \\\e~ nnprat.nr r' 



xT 



is 



the operator exp[z7rA^Jc x | with N a := E xe An xcr and n xa := c xa .c xcr ; the operator c 
replaced similarly, but c x , and cL are unchanged. 

The phonon modes are indexed by j and, for technical simplicity, we assume there 
are finitely many of them, namely v. In some special models, such as the Holstein model, 
defined in (1.9), there is an association between the phonon modes and the lattice sites, but 
such an association is neither required nor assumed. The phonon coordinates are qi, . . . , q v 
and the momenta are pi, . . . ,p u , denoted collectively by q and p. 

The pj's and q/s satisfy the usual commutation relations [qj,Pk] = iSjk an d we shall 
represent these in the usual way as operators on L 2 (R !/ ), the set of square integrable 
functions of v variables, by pj = —id/dqj (with h — 1). 

The most general positive-definite quadratic form can always be put in the normal 
mode form, shown in (1.1), in which the numbers Mj > and Qj > are, respectively, the 
masses and frequencies of the corresponding phonon normal modes. For convenience, we 
explicitly exclude zero frequency modes, which, physically, correspond to center of mass 
translation. The additional potential term V an (q) includes all non-quadratic terms; it is 
completely arbitrary except for the assumption that it is bounded below, i.e., V an (q) > C 
for some number C, and that Y^j=i ^■j^] ( i'j + ^an(q) goes to infinity faster than linearly 
in all directions. 

The real hopping matrix t(q), whose elements are t xy (q), is allowed to be an arbitrary 
measurable function of the phonon coordinates (but not the momenta). An important 
assumption is that t(q) is real and symmetric for each q, i.e., t xy (q) = t yx (q). We also 



assume, for convenience, that Tr |t(q)| is finite for every q. [Here, \t(q)\ = ^t(q) 2 .] We do 
not make any assumption about the relative signs and magnitudes of the hopping matrix 
elements. The reality assumption generically precludes the interaction of the electronic 
orbital motion with magnetic fields. 

The electron-phonon coupling G x (q) is also an arbitrary real function of the phonon 
coordinates that couples the phonons to the total electronic charge at lattice site x. For 
Theorem 1 (existence of singlet ground states) the only assumption about these couplings 
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is lower boundedness of the total phonon potential energy, i.e., we assume that the function 
of q given by 

-Tr |t(q)| - 2 V \G x (q)\ + JVM^ + l-: H1 (q) 






is bounded below. Usually, one assumes that G x (q) is a linear function of q, but we do 
not do so. 

The total spin is a conserved quantity of the Hamiltonian H in Eq. (1.1). The spin 
operators are defined to be the quadratic operators 

S Z '-=\Y,^- n ^ S+ := (S-) f := ]>>W (1-3) 

xeA xeA 

They all commute with the Hamiltonian H. The spin operators satisfy an SU(2) algebra, 
and the total-spin operator is defined to be the corresponding quadratic Casimir operator 
(<Sop) 2 : = {S z ) 2 + t;S + S~ + ^S~S + , with eigenvalues S(S + 1). In particular, we are 
interested in the eigenvalues of the total-spin operator for the ground states of H of the 
electro n-phonon model described in Eq. (1). 

Our main result asserts that the ground state of the model in (1.1) has an S = 
ground state, and that the ground state is often unique. 

THEOREM 1 (Existence of singlet ground states). Assume the previously 
stated conditions on the Hamiltonian H in (1.1). Assume additionally, that the total 
number of electrons, 2N, is even. Then among all of the ground states of H there is at 
least one ground state with total spin S = 0. 

For Theorem 2 (uniqueness of the ground state) additional assumptions are needed. 

(i) The hopping matrix elements t xy are independent of q. We also assume that A is 
connected, i.e., for each x and y in A there are sites x — xq, x±, X2, ■ ■ ■ , x n = y such that 
tx r x i+1 7^ for all < i < n — 1. A bond is said to exist between two sites x and y in A if 
x 7^ y and if t xy ^ 0. 

(ii) All the functions of q appearing in if, i.e., G x (q), V r an (q), are differentiate. [Ac- 
tually, it suffices for them to be locally Holder continuous with densely defined derivatives.] 

(iii) The G x (q)'s are independent. By this we mean that for each point q G R^ the 
v simultaneous equations 

£^U = 0, J = l,2,...,, (1.2) 

xeA ^ J 
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have no common solution other than A x = for all ieA. In other words the (|A| x u) 
matrix dG x /dqj has rank |A| for each q G R^. (Again, it suffices for this to hold only on 
a dense subset of R^.) 

(iv) Every mass, Mj, is finite. 

These conditions hold in many models, e.g., the Holstein model 6 , but conditions (i) 
and (iii) do not hold in the Su-Schrieffer-Heeger (SSH) model 7 . 

THEOREM 2 (Uniqueness of the ground state). //, (i)-(iv) above are satisfied 
then the ground state is unique (and hence a nondegenerate spin singlet). 

Remarks. — (1) Theorem 1 has long been conjectured and is consistent with the intu- 
ition that the exchange of a boson leads to electron-electron pairing. 

(2) The uniqueness theorem establishes that the many-body ground state of H does 
not have any level crossings for a finite system, thereby establishing the result that the 
self-trapping transition from a collection of extended polarons to a collection of localized 
polarons is a smooth crossover, rather than a transition by breaking of analyticity in any 
finite system. 

The proof of these theorems is based upon spin-reflection positivity and is closely 
related to the analogous proof for the Hubbard model already presented by one of us 8 . 
A different proof based upon Perron-Frobenius positivity arguments, was given for one- 
dimensional models 9 , but it does not appear to be readily generalizable to the present 
case. 

The SSH model is the special case of (1.1) in which the hopping matrix elements 
are linear functions of the phonon coordinates and the electron-phonon couplings vanish 
[G x (q) =0]. To be more precise consider the original SSH model on a periodic one- 
dimensional chain 7 

|A| + + 1 |A| 1 

h ssh = ££(*- Wi+i - Qi])(4 + i a c ia + clc l+la ) + - 5>(Q i+ i - Qi) 2 + j^Pf] , 

(1.4) 

with Qi (Pi) the local phonon coordinate (momentum) at site i. Transforming to the 
normal coordinates 



1 |A| 

q ° := ~ 7 fxi ^ Qj 
A| k=i 



cos2nk J -r^ i|A|<j<|A| 

■Uiai ( L5 ) 

sm27Tk J —^ 0<j'<i|A| 
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yields the electro n-phonon Hamiltonian in the form of Eq. (1.1) 

|A| |A|-1 

H ssh = ££[*" T i+M)M +1(J c i(J + clc i+la ) + -J2 l M tflj + mP"] > ( L6 ) 
with 

|A|/2-l |A| 

T,(q)= e (- l )^« in nxr+ E (- 1 )^ cos ^r ' (L7) 

i=o ! ! j=|A|/2 ' ' 

and O 2 = 2k[1 + cos(27Tj/|A|)]/M. Theorem 1 shows that the SSH model always contains 
a spin-singlet ground state, but the ground state is not necessarily unique. We are aware 
of no other rigorous results for the SSH model. 

The Holstein model 6 is a special case of (1.1) where there is one (internal) normal 
mode associated with each lattice site (the index j is identical to the index x) , the hopping 
matrix elements have no phonon coordinate dependence, the electron-phonon coupling is 
linear in the phonon coordinate associated with the lattice site 

Cx(q) = g x q x , (1.8) 

and the potential energy is harmonic [V r an (q) = 0]. The resulting Holstein Hamiltonian is 
HroI = E E **v c Lv + E 9xQx(n x1 + n x[ ) + - ^{M x Sl 2 x ql + j^pI) ■ (1-9) 

a x,yeA xeA x6A x 

The Holstein model has independent couplings if all g x are nonzero. The only rigorous 
result for the Holstein model is that of Lowen 10 for one electron. In this case, it has been 
shown that the ground state is nondegenerate and analytic if the lattice is bipartite. 

In the static limit, where all of the phonon masses become infinite, but the spring 
constant remains finite 

M x -> oo , k x = M x n 2 x = finite , (1.10) 

the phonon kinetic energy terms ^2 X P X /2M X do not contribute to the Hamiltonian (1.9). 
The up- and down-spin electrons become independent and the Holstein model maps onto 
a Falicov-Kimball model 11 with spin one-half conduction electrons and a continuous static 
field q x . Techniques similar to those used in the spinless Falicov-Kimball model 12 may 
be used to show that the ground state is a commensurate charge-density-wave at half 

6 
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filling 13 ' 14 . The static limit of the Holstein model has also been investigated by other 
methods 15 and shown to possess insulating bipolaronic charge- density- wave order at large 
enough coupling. The ground state has also been shown to be nonanalytic in one dimension 15 . 

In the instantaneous limit, where the phonon frequency and electron-phonon coupling 
become infinite, but their ratio remains finite, 

g x -> oo , O x -> oo , -^- = finite , (1-H) 

the Holstein model maps onto an attractive Hubbard model 16 . This mapping is illustrated 
by completing the square in Eq. (1.9) 

H H oi = Yl Yl txyclvCyo ~ g Yl U *( n ^ + n *l) 2 

a x,yEA xEA 

+ ^J2(M x nl[ qx + J ^{n^+n xi }} 2 + ^pl) , (1.12) 

xEA x x x 

with the electron-electron interaction U x defined to be 

U ---A, ■ ai3) 

In this instantaneous limit the remaining electron and phonon terms in the Hamiltonian 
decouple because g x /M x Vt x — > 0, and one is left with an attractive Hubbard model 

H Hub = J2J2 *vLv - E u x n X fn xi + - j2( M Ai + — p d , (i.u) 



X 



cr x,yEA xgA xgA 

with 

vxy • "xy 7) % ®xy • V "/ 

The attractive Hubbard model is already known to have a unique spin-singlet (S = 0) 
ground state for an even number of electrons on a finite lattice 8 . 

In Section II, the proofs of Theorems 1 and 2 are presented for the Hubbard model in 
order to clarify the results of Ref. 8 and to define our notation and methodology. Section III 
contains the proofs of these theorems for the electron-phonon Hamiltonian. A discussion 
of the results follows in Section IV. 

II. ATTRACTIVE HUBBARD MODEL PROOFS 
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Proofs are presented for the results in Ref. 8 in order to clarify the previous work and 
to define the notation and current methodology. We begin with a proof of Theorem 1. The 
attractive Hubbard model Hamiltonian is given by the electronic terms in Eq. (1.14) with 
each U x > and the bar dropped from the hopping matrix (the hopping matrix elements 
have no phonon coordinate dependence here). 

Proof of Theorem 1 for the attractive Hubbard model. — Both the total spin 
operator 5" 2 and the ^-component of spin S z are conserved quantities of the Hubbard model 
(1.14) and restriction can be made to the S z = subspace (without loss of generality), 
because every eigenstate with total-spin S can be rotated into the subspace with S z = 
without changing its energy. Therefore, we can assume there are N electrons of each spin 
(up and down). 

It is convenient to use first quantized notation. We denote the coordinates of the up- 
spin electrons with X which really is an iV-tuple X = (x\, x^-, • • • , xn) where each Xi G A. 
Similarly Y denotes the coordinates of the down-spin electrons. Any wave function \1/ is a 
function of both coordinates \I/(X, Y) and it is antisymmetric in {xi)f =1 and antisymmetric 
in {Vi}iLi- Note that if ^{ x , Y) is an eigenfunction, then so is ^(Y, X) and ^(Y, X)*. (It 
is here that the condition that the hopping matrix elements t xy are real is used.) Instead of 
considering an eigenfunction \I/(X, 1") it is convenient to consider ^(X, Y) + ^>(Y, X)* and 
i[ty(X, Y) — \I/(Y, X)*]. In other words we can, without loss of generality, assume that our 
eigenfunction, viewed as a matrix indexed by X and Y, is self adjoint (but not necessarily 
real) . 

-$(X,Y) = *(y,X)*. (2.1) 

The dimension, d, of this matrix \1/ is 

Any self-adjoint matrix can be expanded in terms of its eigenfunctions. Thus, 

d 
tf (X, Y) = J2 w a MX)MY)* (2.3) 

a=l 
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where the </> a 's are an orthonormal set of functions (but antisymmetric in their argument 
X), and the w^s are real numbers. Our aim is to show that the w a 's can all be chosen to 
be nonnegative. This will conclude the proof because it implies that 

J2^(X,X)>0, (2.4) 

x 

which implies that tfj(Xo,Xo) is positive for at least one Xq. This means that the wave 
function does not vanish when the up-spin electrons and the down-spin electrons are at 
precisely the same locations - - and this is a singlet state. Thus i/j has a nonvanishing 
component in the S = sector. 

To show that the w^s can be taken nonnegative, let us write out the energy using 
the decomposition of \l/ in (2.3). One easily computes 

d |A| d d 

{^\H\^) = 2^wl((l> a \K\(l> a ) -^Ux^^WaW^alL^)^ (2.5) 

a=l x=l a=l /3=1 

d 

(*|*> = X>2<0«|0«> (2.6) 

a=l 

where the <i-dimensional matrices K and L x are defined as follows: K is the first-quantized 
version of ^2 xy t xy c x c y (no spin here) and L x is the first-quantized version of n x . More 
explicitly, the matrix elements appearing in (2.5) and (2.6) are constructed in the following 
manner: The inner product between two arbitrary vectors 4>i{X) and 4>2{X) is 

(0 2 |0i) = J>2pO*0i(*)- (2-7) 

X 

The kinetic energy matrix elements satisfy 

JV 

(0 2 |-K10i) = X^XIX^ 2 ^ 1 '---'^' ■■■iX N )*(j) 1 (x 1 ,...,y, ...,x N ) (2.8) 

xy j=l X 

where the argument of 02 agrees with the argument of 4>i everywhere except at the j 
index where the site index is x for <p2 and y for (pi . The number operator matrix elements 

satisfy 

TV 

(<h\L x \<l>i) = Y,MX)*MX)Y, S *i'* ( 2 - 9 ) 

X j=l 

9 
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For the purpose of Theorem 1, the explicit values of these matrix elements are unim- 
portant. The only thing one has to note about (2.5) and (2.6) is that replacing every w a by 
\w a \ cannot increase the energy. The first term in (2.5) and the inner product of (2.6) stay 
the same, while the second term in (2.5) can only improve (if it changes at all). Thus if \I/, 
given by (2.3) is a ground state of the Hubbard model then so is |\I/| which is constructed 
by replacing w a — > \w a \ in (2.3). Note that I^K-X", Y) is not equal to |vP(-X", Y)\ in general, 
but corresponds to |\I/| = yW^ in the sense of matrices. Q.E.D. 

In order to prove Theorem 2 we first need to establish a lemma and in order to state 
the lemma a definition is needed. If X = (xi, . . . , xn) we define the operator II X to be 

U x :=L X1 L X2 ...L XN (2.10) 

Since the different L x s commute, the ordering of factors in (2.10) is unimportant. In second 
quantized notation U x = n Xl . . .n XN which shows that Ii x is an orthogonal projector, i.e., 
U x = (II x )t and (II X ) 2 = U x . It is also obvious that II X is a one- dimensional projector. 
Furthermore, if X and X' differ only by a permutation then Ii x = U x . The matrix 
elements of U x [analogous to (2.9)] satisfy 

(0 2 |n x |0i) = MxyMx)N\. (2.ii) 

[The AH in (2.11) may appear mysterious, but it is not. The reason is that if <p\ and fo are 
normalized vectors concentrated at X and all of its AH permutations and at no other X, 
then \(pi(X)\ 2 = \4>2{X)\ 2 = 1/Nl. Thus, both sides of Eq.(2.11) are equal in magnitude 
to 1.] 

LEMMA (Connectivity of the single-spin configuration space). Assume 
that the lattice A is connected as explained above. Then the single-spin configuration 
space is connected by the kinetic energy matrix K. That is to say, given points X = 
(xi, . . . , xn), Y = (yi, . . . , y^v) in the single-spin configuration space, there exists a chain 
of m points [Y = X m ,X rn _i, . . . ,X%,Xi = X} in the configuration space, such that the 
product of matrix elements satisfies 

u x mKIl x m _ 1 u X3 KU X2 KU Xl ^ 0, (2.12) 

withX 1 ^X 2 ^X 3 ^...^X m . 

10 
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Proof: We first consider a geometric question. Place iV unlabelled markers on the 
points (xi, . . . , xn) of A. The goal is to move these markers, one at a time, across bonds of 
the lattice A to a final set of points (j/i, . . . , x/n) in such a way that at no step is there ever 
a doubly occupied site of the lattice. We want to emphasize that it is not necessary that 
the marker that was first at x\ ends up at y\, we require only that in the final configuration 
the sites (j/i, . . . , y^) are occupied by some marker. 

To do this we apply the following algorithm repeatedly — at most A" times. Look for 
the smallest i such that the site yi is unoccupied. Look for the smallest j such the site Xj is 
not in the set {y\, . . . y/v}. We will move a marker from the point Xj and establish a marker 
at the point yi in such a way that the other occupied sites of the final state are identical 
to the other occupied sites of the initial state. (Once again, we point out that the marker 
originally at Xj need not end up at yi, and the other markers may be moved in this process.) 
To achieve this we choose a connected path P in the lattice A, P — (xj = z\, z%, . . . Zk = yi) 
from site Xj to yi. Such a path exists by hypothesis. If there are no markers on the sites 
Z2, ■ ■ ■ , Zk-i, then we simply move the marker at Xj along the path P to y^. Suppose on the 
contrary, that there are some other markers on this path P. Let / be the largest number 
such that z\ has a marker on it. Then simply move this marker along the path to the 
site yi, thereby achieving two things: a marker on yi and one less marker along the path 
P. We then move in turn each marker on the path P to the location of the previously 
moved marker. This completes the description of the algorithm and answers the geometric 
question. 

To prove the lemma itself, we first note that (4>2\n Xl KIl X2 \(j)i) = for all X\ and 
X2 unless there is a permutation of X\ such that after the permutation there exists some 
1 < J ' < A" such that x\ = x\ for all i 7^ j. In this case 

(^lU^KU^l^) = t ab MXirMX 2 )Nl (2.13) 

where x] = a and xj = b. 

The configurations Xi,X 2 , ■ ■ -X m used in Eq. (2.12) will be the X's determined by 
the sequence of moves in the geometric discussion above. (Note that although the markers 
were indistinguishable there we can, if we wish, put numbers on them. In general the 
final state in this case will not be identically the state Y = (yi, . . .yn) but will be some 



11 
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permutation Y of the final state. On the other hand H Y = H Y because the L x s commute.) 
What remains to be shown is that the operator 

U Xm KU Xm - 1 '--U X2 KU Xl = (IL Xm KIL Xm - 1 )(n x — x KU Xm ~ 2 ) ■ ■ ■ (U X2 KU Xl ) (2.14) 

is nonvanishing. Let 0i, . . . , 4>d be the orthonormal basis proportional to antisymmetrized 
delta-functions in the configuration space. From (2.13) we see that when we expand (2.14) 
in these intermediate states, that a factor such as ((f>i\Il Xk KIl Xk - 1 \(f)j) can be nonzero for 
only one i and j, and this contribution is equal to t a h (up to an overall sign) where a and 
b denote the two indices where Xk and Xk-i differ. So, in short, there will be exactly one 
matrix element of ((j)i\H Xm KH Xm - 1 . . . H X2 KH Xl \<f>j), and it will be a product of t X yS all 
of which are nonzero. Q.E.D. 

We turn now to the proof of Theorem 2 for the Hubbard model. The four assumptions 
(i)-(iv) for the electron-phonon model simplify to the two assumptions: (a) A is connected; 
(b) every U x is positive (U x > 0). 

Proof of Theorem 2 for the Hubbard model. Suppose \&i and \&2 are orthogonal 
ground states. We can assume both \&i and \&2 are Hermitian. Then \I/ = \I/i + A\&2 is a 
ground state for all real A. Moreover \1/, if viewed as a matrix, cannot be either positive 
or negative semidefinite for all real A. Hence for some choice of A we have that the two 
ground states \l/± := |(|\1/| ±\1/), are both nonzero, and are both ground states, since |\I/| is 
a ground state from Theorem 1. In particular \1/ + satisfies the matrix Schrodinger equation 

|A| 

K^ + + ^ + K-J2u x L x ^+L x = e^+, (2.15) 

x=l 

and is a positive semidefinite matrix. 

We define 7i_|_ to be the range of the matrix \1/_|_, which is a subspace of H = C d . We 
also define H± to be the orthogonal component of H+ and n^ to be the projector onto 
7i±. By assumption H.+ and 7i± are both nontrivial. 

Multiply the Schrodinger equation (2.15) on the left and on the right by n^ to yield 
Y, x ^ r xH_ L L :r \l/ + L x n_ L = since H±_^ + = ^+T1± = 0. This implies further that 

tf+LJI ± =0 (2.16) 



12 



JKFEHL - 29/June/94 

since every U x is positive and \I/-|- is a positive semidefinite matrix. Eq. (2.16) implies that 
the matrix L x does not connect the subspaces H+ and Hj_. Now multiply (2.15) on the 
right by II j_ and use (2.16) to discover that 

*+Mi = 0. (2.17) 

Both (2.16) and (2.17) show that 7i+ and 7i±_ are invariant subspaces of the matrices K 
and L x . 

We will now use the lemma to show that 7ij_ is trivial, thereby establishing a con- 
tradiction. The operators U x defined in (2.10) are one- dimensional projectors in TC and 

satisfy 

(A^!)- 1 ^n x = l. (2.18) 

x 

Thus V + = U x ty + is nonzero for some X. Likewise, if E 7i± then V 1 - = U Y <fi ^ for 

some Y. By the lemma, the rank one operator in (2.12), call it A, is nonzero and satisfies 

(V- L ,AV + ) 7^ 0. This is a contradiction since A, being a product of L^'s and K, has H+ 

as an invariant subspace. Q.E.D. 

III. ELECTRON-PHONON MODEL PROOFS 

The proof of Theorem 1 employs the same methods as utilized for the Hubbard model 
in the preceding section, but now all of the expansion coefficients and basis functions have 
an implicit q dependence. The variational principle that shows that the ground state 
includes a spin-singlet state will arise from the kinetic energy terms for the phonons, since 
there is no direct electron-electron interaction here (however, see the remark at the end of 
the proof of Theorem 2 about including attractive electron-electron interactions) . 

Proof of Theorem 1 for the electron-phonon model. Any wavefunction \I/(q) 
can be thought of as a matrix-valued function of q G R^. By 577(2) invariance of (1.1) 
we can, as in the preceeding proof, assume that there are equal numbers of up-spin and 
down-spin electrons, namely N of each kind. The dimension d of the matrix \I/(q) is then 

As before, by taking transposes and complex conjugates, we can restrict our discussion to 
the case where \I/(q) is a Hermitian matrix for all q (it is here that we use the condition 
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that the hopping matrix elements t xy are real because, without this condition, the complex 
conjugate of \1/ will generally have a different energy from that of \I/). Note that ^(q) is 
Hermitian, but it need not be real. 

We can write the Schrodinger equation for \& in the following generic way 

v 1 rfl 

This equation is to be understood in the following sense: \I/(q) is a matrix- valued function 
and so are its second derivatives; V(q) is also a matrix- valued function which is self-adjoint 
for every q; the two terms V\E'-|-\E'V A (with matrix multiplication being understood) include 

V 

all of the terms in 7i besides the phonon kinetic energy, Yl P^j/^Mj; of course e is the energy 

eigenvalue. 

Associated with (3.2) is a variational expression which we can write as £(\I/)/(\I/|\I/). 
The denominator has the form 



(#|#) = / Tr^(q) / dq. (3.3) 

The numerator is 

*w = f n Aip 2^ Tr[ ^* (q)]2 + 2M(q)$(q)2 1 dq ' (3 - 4) 

where dj denotes the partial derivative d/dqj. 

Our strategy is to replace the matrix ^(q), for every q, by its absolute value in the 
matrix sense, i.e., 



|*(q)| = V*(q) 2 - (3-5) 

We note that the norm satisfies (^l^) = (|\I/| |\I/|), and that the V^ 2 term is evidently 
unchanged. In the appendix we prove that this substitution does not increase the integral 
of Tr[<9j\I/(q)] 2 . There are some nontrivial technical issues here caused by the fact that 
dj^>(q) and <9j|\I/(q)| may only be distributional derivatives, but these issues are fully 
resolved in the appendix. 

Since, by definition of the ground-state energy, £(^) cannot decrease when \1/ is re- 
placed by |\I/|, we conclude that |\1/| is also a ground state of the Hamiltonian. As in the 
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Hubbard model proof, we conclude that |^/(q)| has a nonzero projection onto the S = 
subspace, i.e., if every diagonal element of the matrix |^/(q)| vanished for almost every 
value of q, |^/(q)| would be the zero function, which it is not. Q.E.D. 

Proof of Theorem 2 for the electron-phonon model. Suppose there are two 
ground states ^i(q) and ^(q) which are Hermitian and linearly independent. Then 
\l/i + X^2 is a nonzero ground state for every real A, and as A is varied from — oo to +00 
there will be values of A for which \& = \&i + A\&2 has both a negative and positive spectrum 
for a set of q's of positive measure. Then |\I/| is also a ground state, as are \I/+ := -|(|\I/|+\I/) 
and \l/_ := |(|\I/| — \&). Note that W, |W|, \1/_|_ and \I/_ are all nonzero functions that satisfy 
the Schrodinger equation (3.2). Indeed \1/ = \I/ + — \I/_. Moreover, both \I/_|_(q) and \I/_(q) 
are positive semidefinite matrices for all values of q. From the fact that the matrix valued 
function V(q) appearing in (3.2) is differentiable, elliptic regularity theory applied to the 
Schrodinger equation tells us that \I/(q) is twice continuously differentiable [actually only 
Holder continuity of V(q) suffices for this conclusion]. 

At each point q, the vector space C d , on which \l/(q) operates, is naturally the direct 
sum of three subspaces (some of which might be empty). These subspaces are denoted 
by 7i + (q),7i_(q), and 7io(q)- H + (q) is the spectral subspace of \& + (q), i.e. consists 
of all linear combinations of the nonzero eigenvectors of \I/_|_(q); 7i_(q) is the spectral 
subspace of \&_(q); and 7io(q) is the orthogonal complement of 7i+(q) © 7Y_(q), i.e., the 
subspace of all linear combinations of the zero eigenvectors of \l/_|_(q) and \l/_(q). Of course 

C d = ft + (q)©ft (q)©H-(q). 

Let <io(q) denote the dimension of Ti.o(q) and let do denote the minimum of {do(q) : 
q G R^}. Since there are only finitely many values for <io(q) there is a point qo G R^ for 
which do(qo) — do- By definition the matrix ^(qo) has do zero eigenvalues, and the positive 
eigenvalues are separated from zero by a gap A. Since \1/, and hence the eigenvalues of \1/, 
are continuous, there is some ball B' centered at qo with radius r such that the dimension 
of do(q) < do for every q G B' . Since, however do is the minimum of Gfo(q), we conclude 
that do(q) = ^0 for every q G B' . 

Now let us study the contour integral 

If 1 

n ± (q) := / -— dz (3.6) 

±KHJ 27riy r *(q)-z l ' 
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where the contour F runs from — oo just below the negative real axis, goes vertically upward 
when the real part of z = A/2 and returns to — oo just above the negative real axis. 

First we observe that n_i_(qo) is the projector onto the orthogonal complement of 
7i_|_(qo), namely onto Ti.±(qo) = ^(qo) © ^-(Qo)- Furthermore, the zero eigenvalues of 
\I/_i_(q) do not move from zero as long as q G £>', as we have just seen above. Therefore, 
n_i_(q) continues to be the projector onto Ti.±(q) as long as the positive eigenvalues of \l/(q) 
are greater than A/2, and hence do not intersect the contour F. Since \1/ is continuous we 
conclude there is a smaller ball B C B' centered at q in which n_i_(q) continues to be the 
projector onto Hj_(q). Since \I/(q) is twice continuously differentiate, it is a trivial matter 
to show that we can differentiate under the integral sign in (3.6) and conclude that n_|_(q) 
is a twice continuously differentiable matrix-valued function on B. 

Now we compute some derivatives in this ball B. We start with the observation that 
\E f _l_(q)ilj^(q) = for all q in R,". The following identities hold in B: (we suppress the q 
dependence) : 

(^■*+)n ± + * + ^n ± = (3.7) 

(d^+)u ± + 2(d j ^+)d j u ± + # + dJn ± = o. (3.8) 

If the Schrodinger equation (3.2) for \I/_|_(q) is multiplied on the left and on the right by 
n_L(q) we have (since ^ + n_L=0 = n ± * + ) 

E^n ± (dJvi> + )n ± = o, (3.9) 

or, combining (3.9) with (3.8) we discover that 

J2^r u ^ d ^+^ U ^ = °- ( 3 - 10 ) 

3=1 J 

Now multiply (3.7) on the left by djFi±/2Mj and sum over j. The first term vanishes 
because of the adjoint of (3.10). The second yields 



E 9XT^ n -L*+ a i n -L = °- t 3 - 11 ) 



2M. 

3 = 1 3 
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Since \l/_|_(q) is a positive semidefinite matrix (and all Mj < oo), we conclude from (3.11) 
that 

^+djU ± = 0, j = l,...,u (3.12) 



for all q G B. Eq. (3.12) states that the range of the derivatives of II j_ lies in 7i±, or 



diU± = U ± d j U ± 



J^-D 



(3.13) 



and therefore that each term in (3.7) separately vanishes, i.e., 



(d^ + )Il ± = 0. 



(3.14) 



Differentiating (3.14) yields 



(0?tt + )n_L + (<9^+)<9,n ± = o. 



(3.15) 



But (djV + )djIL± = (dj^!+)U±djU± = from (3.13) and (3.14), so we finally conclude 
that the range of the second derivative of \1/_|_ lies in 7Y+, or 



(d^+)U ± = 



(3.16) 



for all q G B. 

If we multiply the Schrodinger equation (3.2) for W_|_ on the right by II j_ and use the 
identity (3.16) we find 



*+(q) 



|A| 



K(q) + J2L x G x (q) 



x-l 



n ± (q) = o. 



(3.17) 



This says that the matrix W(q) := K(q) + Yl x =i L x G x (q) does not connect 7i_|_ with H±. 
Differentiating (3.17) shows that 



(djV + )Wn ± + ^+(d J W)U 1 _ + ^ + WdjU± = 0. 



(3.18) 



But we know that the range of the derivatives of \1/+ lies in 7i+ and that of the derivatives of 
U± lies in 7ij_, so the first and last terms of (3.18) vanish because W does not connect 7i+ 
with 7i± . Since the hopping matrix elements do not have any q dependence by assumption 
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and the derivative of the electron- phonon couplings is a rank |A| matrix by assumption, 
we find 

^+L x U ± =0 x=l,...,|A| (3.19) 

and, from (3.17), 

^ + KU ± = 0. (3.20) 

These two identities imply that both 7i + (q) and 7Yj_ (q) are invariant subspaces of the 
matrices K and L x for all q G B. Exactly as in the Hubbard model proof, we conclude that 
for every q G £> one of the two alternatives holds: either H+(q) = {0} or 7Y_(q) = {0}. 
Since the functions *&+ and \l/_ are continuous, the set on which \I/_|_ is nonzero is open 
and the set on which \l/_ is nonzero is open. Therefore B contains an open set in which 
either \I/_|_(q) is identically zero or \l/_(q) is identically zero. However this cannot happen 
since \1/_|_ and \l/_ are eigenfunctions with locally bounded potentials V^(q) and therefore 
satisfy a unique continuation theorem 17 . That is to say if \I/-|- vanishes in some open set, it 
vanishes in all of R" contrary to our original assumption that \I/-|- is not identically zero. 
Q.E.D. 

Remarks. — (1) Both Theorems 1 and 2 continue to be valid if any attractive Hubbard 
model terms are added to the Hamiltonian in (1.1). In other words, the electron- phonon 
model and Hubbard model results do not interfere with each other as long as they are 
attractive. 

(2) The SSH model is shown here to have at least one spin-singlet state among its 
ground states. The ground state is not shown to be unique here. The difficulty for the 
uniqueness proof enters in the derivative of the matrix W in Eq. (3.18). 

IV. DISCUSSION 

The results presented in this contribution hold only for an even number of electrons 
in a finite system. In this case, the only nonanalyticities that can enter in properties of 
the ground state occur when there is a ground-state level crossing. In the cases where the 
ground state can be shown to be unique, there are never any level crossings, so that any 
"transition" of the electron-phonon ground state from a collection of delocalized polarons 
to a collection of self-trapped polarons is not a sharp transition, but is a smooth crossover. 
We are not prepared to prove any statements about the thermodynamic limit here. 
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One can ask if these results will survive if a magnetic field is turned on. The answer 
in general is no because the variational argument presented in Theorem 1 no longer holds 
if there are interactions with the local electronic spin [as one would have if one added a 
Zeeman coupling to the Hamiltonian in Eq. (1.1)]. One can also investigate the effect 
of a magnetic field on the electronic kinetic energy. The hopping matrix elements are 
always assumed to be real, and therefore can only incorporate flux phases that correspond 
to entrapped fluxes that are integral multiples of n. If the hopping matrix elements be- 
come complex, the entire methodology incorporated here fails, and there are no rigorous 
statements that we can make about this case. 

In conclusion, we have presented a proof that the ground state of a general class of 
electro n-phonon Hamiltonians must include a state that is a spin singlet. The phonons 
can interact with the electrons in two different ways — the phonons interact with the local 
electronic charge and the phonons modulate the electronic hopping integrals. The phonon 
coordinate dependence of the both electron-phonon couplings and the hopping matrix ele- 
ments is arbitrary. The phonons can be optical modes or acoustical modes, and can contain 
anharmonic couplings. The hopping matrix is always assumed to be real and symmetric. 
The Su-Schrieffer-Heeger model, the Holstein model, and the Hubbard model all fall into 
this general class. In the case where the hopping matrix contains no phonon coordinate 
dependence, the lattice is connected, the electron-phonon couplings are independent, and 
the inverse phonon masses are all positive, the ground state has also been shown to be 
unique. 
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APPENDIX: ABSOLUTE VALUE DECREASES KINETIC ENERGY 

We shall prove here that replacing a matrix- valued function, \I/(q) by its absolute 
value, |^/(q)| = ^/\I/t(q)\I/(q), decreases each component of the kinetic energy 



TjW := J IV[^*t( q )^*(q)]dq, (A.l) 



R'- 



when \1/ is Hermitian. (Here dj = d/dqj and f denotes adjoint.) 

Before going into the technicalities, let us give a heuristic discussion to motivate the 
truth of our assertion. If we write, in Dirac notation, 



*(q) = X> a (q)|/ a (q))</ (q)|, ( A2 ) 



a=l 



where d is the dimension of the matrix (in our proof later we shall generalize to d = oo), 
the functions w a are the real eigenvalues and the |/ Q (q)) are the q-dependent orthonormal 
eigenfunctions of the Hermitian matrix ^(q). Supposing everything to be differentiable 
we can compute 

d 
0,-tf(q) = ^^ a (q)|/ a (q))(/ Q (q)|+«; Q (q)|/ a (q))^ Q (q)|+w a (q)b a (q))(/ Q (q)|, (A3) 

a=l 

where \g a (q)) = dj\f a (q)). Since (/ a (q)|/g(q)) = $ a /3, we have that (/ Q (q)|^(q)) + 
(fl'a(q)l//9(q)) = 0- Thus, if we square (A. 3) and take the trace we have 

d d d 

a.-\ 01=1/3=1 

(AA) 
From this we see that replacing w a (q) by |tu Q (q)| can only decrease the last term on the 
right. The second term does not change. The first term also does not change since, for 
any real, differentiable function w(q), it is a fact that [djw(q}] 2 = [dj\w(q)\] 2 in the sense 
of distributions. 

Notice that this heuristic discussion gives a pointwise inequality Tr[<9j\I/(q)] 2 > Tr[<9j |\P(q) 
In our rigorous discussion we shall content ourselves with an inequality for the integral 
(A.l) — which is sufficient for our purposes in this paper. 
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We begin the rigorous discussion with some definitions. Let 7i be a fixed, separable 
Hilbert space (for the purposes of our paper H is finite dimensional, but there is no need 
for this restriction here). Let B denote the Hilbert-Schmidt operators on 7i and, for A G B, 
let 

||A|| := {TrA^A} 1 / 2 

denote its Hilbert-Schmidt norm. We also define 



\A\:=VaTA and \A^\ := J~AA^ 

and note that |||A||| = |||At||| = \\A\\. 
A map 

* : R^ -> B 

is said to be measurable if the function /a(q) : = ||^(q) ~~ ^-11 is (Lebesgue) measurable for 
every Hilbert-Schmidt A. It is not hard to prove that \& is measurable if and only if every 
matrix element (v,^(q)w) is a measurable function for every fixed v and w in 7i. \1/ is 
said to be in L 2 (R^; B) if \& is measurable and if / G L 2 (R^), i.e., j fo(q) 2 < oo. 

The map \1/ is said to be in i7 1 (R^;i3) if \& G L 2 (R^;£>) and if there are maps 
<9i\I/, 82^, • • • , d u ty for which the following holds 
(i) Each djV : R" -»■ B is a map in L 2 (R^; B) 
(ii) For each infinitely differentiable map, : R" — ► £>, of compact support with derivatives 

dj(j), we have the relation 

r TV[a-0(q)*(q)]dq = - / Tr[0(q)&-*(q)]dq (A5) 



R" R^ 



for each j. (Note: to say that <j) is differentiable means that for every q G R", lim e ^o ||s X [</ ) (q+ 



ee_j) — 0(q)] — <9j0(q)|| = with e^ being the unit vector in the j direction. Since <f> has 
compact support this limit is uniform in q.) Clearly 8j^ = (dj^y. 

It is easy to verify that H 1 (R 1 '; B) is a Hilbert space with inner product 

(*.*^/^W(q) + X>*W'( q )]d<, (-4.6) 



R' 
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Clearly, \1> is in H 1 (or in L 2 ) if and only if ^ is in H 1 (or in L 2 ). This class, i/ 1 (R !/ ; i3), is 
precisely the class needed for quantum mechanics, i.e., so that the variational energy Tj(^) 
can be defined and so that the norm J Tr\I/t\I/ can be defined. For matrix-valued functions 
(i.e., 7i is finite dimensional) the properties of measurability, being in L 2 and being in 
H 1 are just the ordinary meaning of these properties applied to each matrix element of \1/ 
considered as a function on R^. 

We are indebted to Jan Philip Solovej for very considerable help with the following. 

THEOREM. Let ^ be in ^(R^S). Then \V\ and |tft| are i n H X {W;B) and, for 
each j 

2T j (*)>T J m)+T j (\&\). (A.7) 

In particular, if \I/(q) is self-adjoint for all q then 

T^yT.m). (A8) 



Remark: If \1/ is not self-adjoint it is quite possible that Tj(^>) < Tj(|\l/|). But then 
T i (*)>T i (|*t|). 

Proof: It suffices to prove the theorem when \I/(q) is self-adjoint for all q. To see 
this, consider the Hilbert space H2 = H. ® Ti and replace \1/ by the self adjoint operator 

* 2= (£ * t )- Thfin l*2l=( l o' |*t|) and (A.7) becomes ^(^^T.d^l). 

The first task is to show that |\1/| is measurable. Our definition of measurability given 
above is the standard one, namely the inverse image of an open ball in the Banach space 
B is measurable. Now the map \1/ — > |\I/| is continuous in the B norm and hence, by a 
standard result, |\1/| is measurable. 

Let vi,V2-.., be an orthonormal basis for 7i. Matrix elements (v a ,Avp) will be 

denoted by A a p. It is easy to verify that dj{4> a p) = (dj<p) a/ 3 in the classical sense if 

4> is differentiable. If \1/ has distributional derivatives then dj($ a p) = (djtyotp in the 

ordinary distributional sense. [Simply take 0(q) = h(c\)\vp)(v a \ in (A. 5), with h being 

an ordinary infinitely differentiable function of compact support.] Another important 

preliminary remark is that by Fubini's theorem J^ J RI , K ot /3(q)dq = J R „ ^2 K a p(q)dq for 

a, P a,/3 

any nonnegative, measurable functions K a p(q). 
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For any \1/ in H 1 (R"; B) (including 0, as a special case) we define <9 J£ \I/(q) = e 1 [^>(q+ 
sej) — \I/(q)]. By the fundamental theorem of calculus for distributions 

[d je *(q)] afi = d je Mq) a p] = [ d j [*(q + tee j ) al3 ]dt= f [d^{q + Uej)] a pdt (A9) 

Jo Jo 

Thus, by Fubini's theorem and Schwarz's inequality, 

/ Tr[c» J£ vl/(q)] 2 dq = V f ( f [d^(q + tee^Udt) dq 
Jn v p Jr v Uo ) 

<T [ I [d j V(q + tee j )]i f3 dtdq= f Tr[^*(q)] 2 dq. (A10) 

Q/3 Jw Jo Jr.* 

Therefore, dj e ^> is uniformly (in e) bounded in L 2 (R !/ ; B). 

Now we are ready to study |\I/|, and we begin with a little lemma. If N and M are 
self-adjoint Hilbert-Schmidt operators then 

Tr(N-M) 2 >Tr(\N\-\M\) 2 . (All) 

This is equivalent to TriVM < Tr\N\\M\. If we write 2N± = \N\ ± N and 2M± = 
\M\ ± M, we have that N± and M± are positive semidefinite, and our requirement now 
reads Tr [M + iV_ +M_N + ] > 0. But this is true because TrM + iV_ = TrM* /2 iV_M* /2 > 0, 
etc. 

From (A. 11) we have that Tr[<9j £ |\I/|] 2 < Tr[dj £ ^] 2 , from which we deduce that dj- e |\I>| 
is uniformly bounded in L 2 (R U ;B). More precisely, by (A. 10) and (A. 11), 

f [a je |vl>|(q)] 2 dq< / Tr[a,vl>(q)] 2 dq. (A12) 

Since L 2 (R^; B) is a Hilbert space, the boundedness shown above plus the Banach-Alaoglu 
theorem implies that there is a sequence £i,£2, • • • , tending to zero such that 

fyeM^Pj, (A13) 

as n — > oo, where Pj is a map in L 2 (R^;£>) and where the convergence is in the weak 
sense. Since, dj S (j)(q) converges to dj(p(q) uniformly (in Hilbert-Schmidt norm) we can 
write (A. 5), using (A. 13), as 

/ Tr[0Pj] = lim / Tr[0^ e |*|] = - lim / Tr[<9j £ 0|#|] = - / Tr[<9,-0|*|] (A14) 

JR- ^"JR" ' e ^°jR^ JR" 
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The first equality in (A. 14) is the weak convergence of dj e ^>; the second is just a trivial 
change of variables in the (/-integration; the third is the uniform convergence of dj e <p. 
Equation (A. 14) holds for every (p. By uniqueness of the distributional derivative for 
ordinary functions [and choosing 0(q) = h(q)\vp)(v a \ as before] we conclude that (Pj) a /3 = 
(<9j|\I/|) Q/ 3, and hence Pj — <9j|W|. However, norms cannot increase under weak limits and 
thus (A.8) follows from (A.12). Q.E.D. 
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